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We review the possibility to use the renormalons emerging in the perturbation series of the 
1^ twist- 2 part of the nonsinglet structure functions Fi, F%, F3, and g\ to make approximate 

Q> ■ predictions for the magnitude of the appertaining twist-4 corrections. 
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u 1 Introduction 

The precision of deep inelastic scattering experiments nowadays allows for the disentanglement of 
genuine twist-4 corrections to unpolarized structure functions. While the twist-2 parton density 
has a simple probabilistic interpretation in the parton model, genuine twist-4 corrections can 
be seen as multiple particle correlations between quarks and gluons. For example the twist-4 
correction to the Bjorken sum rule f^ can be interpreted as the interaction of the induced 
color electric E a A and color magnetic fields B a A of a single quark and the corresponding nucleon 
remnant [IJ: 

( P S\-B A j A + A xE A y\pS) = 2m 2 N fWS° , (1) 

where is the nucleon mass, p the momentum and S the spin of the incoming nucleon. Despite 
of the capability to isolate at least for the unpolarized case the twist-4 contributions to structure 
functions in experiments , on the theoretical side the description of higher- twist contributions 
is still not satisfactory. Even though the operator product expansion is possible for inclusive 
quantities its technical applicability is restricted to the lowest moments of the structure functions 
. The best and physically most stringent approach to estimate the magnitude of the matrix 
elements of the twist-4 operators, the lattice gauge theory, suffers from the problem of operator 
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mixing which still has to be solved. In this approach up to now higher-twist corrections were 
calculated only in cases where such a mixing is forbidden by some quantum numbers, e.g. as in 
the case of the twist-3 correction to the Bjorken sum rule Qj. In total, we are far away from get- 
ting hold of the twist-4 contribution to all moments and consequently from a description which 
gives the x-dependent contribution of the twist-4 correction to the whole structure function. 

The renormalon approach cannot really cure this problem because of its approximative and 
partially very speculative nature, but it offers a well defined model. This model was applied to 
the nonsinglet structure functions either in terms of the massive gluon scheme or the running 
coupling scheme which in general have an one to one correspondence ||. Empirically, it has 
proven to be quite successful to reproduce at least in the large-x range the x dependence of the 
twist-4 correction to the unpolarized structure functions |6|, [7|, |8| . The success of this model is 
nowadays referred to as the 'ultraviolet dominance' of twist-4 operators M. 

In the following first two sections we will shortly outline the definition and the basic features of 
the renormalon in QCD and the technique of naive nonabelianization (NNA), which allows to 
reduce the calculational effort necessary to the mere summation of vacuum polarization bubbles. 
We also comment on the very important question of the scheme dependence of the renormalon 
and in the last section we want to discuss the renormalon contribution to the structure functions 
gi, Fl, F 2 , and F 3 in detail and compare it to the measured data. For a review and references 



to the classical papers on renormalons see [10]. 



2 What are Renormalons? 



The operator product expansion predicts that a moment of an unpolarized structure function F 
(F = i*2, Fl, Fs) admits up to 0(1/Q A ) the following general form: 



M n (F) := f dxx n F{ 
Jo 



x 



f~i twist=2 i 



M 2 C' n b^ ist=2 + C^b^ st=4 + C n 2) b^ st=A + .. 



(2) 



Each Wilson Coefficient (e.g. C n ) can be written as an asymptotic series in the strong coupling 
constant. 
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C n = C n (Q 2 /^ 2 , a s ) = J2 C m, n aT + AR; a s = -± . 

m=0 47F 



mo corresponds to the minimal term of the series and is given by the condition 
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AR = C n ,m a 
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(3) 



(4) 



In a certain approximation the Wilson coefficient can be evaluated to all orders, i.e. when we take 
the numbers of fermions to infinity, iVp — > oo. One then can derive a closed formula for C mj „ for 
all m and n. Only the flavor nonsinglet case has been considered so far. The Borel-transformed 
series 
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with p = —j^^Nf (Np —> oo), has the pole structure: 

/ n 2 p -C\ s ( yJRz n IRi n UVi n UV 2 \ 

B[c„](s) = + £_ + ik_ + J?=_ + £_ + ...) , (6) 

reflecting the fact that the QCD-perturbation series is not Borel summable. The pi 3 and 
Pn Rj define residua of ultraviolet- and infrared-renormalon poles, respectively. The existence of 
infrared-renormalon poles makes the unambiguous reconstruction of the summed series from its 
Borel representation impossible. C is the finite part of the fermion loop insertion into the gluon 
propagator equal to —5/3 in the MS scheme. The Borel representation B[C n ] can be used as the 
generating function for the fixed order coefficients 



d m 

Cn,m = (3^B[C n }(8) 



(7) 

s=0 



which implies for large m the following asymptotic behavior of the coefficient functions 



/ 2 -C\- 2 | 
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The infrared renormalon pole nearest to the origin of the Borel plane dominates the asymptotic 
expansion. Its name originates from the fact that in asymptotically free theories the origin of 
the factorial divergence of the series can be traced to the integration over the low momentum 
region of Feynmann diagrams: 

with 2 

a(fc)= l + /?o«(Q 2 )ln(P/Q 2 ) • (10) 

The uncertainty of the asymptotic perturbation series can be either estimated by calculating the 
minimal term in the expansion (|3|) or by taking the imaginary part (divided by 7r) of B[C n ](s) 

Q: 



1 1 f°° ( —s \ 

AR = -I- dsexp[-—)B[R](s) 

7T fj Jo \Poa s J 

The undetermined sign is due to the two possible contour integrations, below or above the pole. 
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In a physical quantity the infrared-renormalon ambiguities have to cancel, against similar 
ambiguities in the definition of twist-4 operators such that in principle only the sum of both 
is well defined. Operators of higher-twist and therefore higher dimension may exhibit power- 
like UV divergences, and the corresponding ambiguity can be extracted as the quadratic UV 
divergence of the twist-4 operator. Because of this cancellation, the IR renormalon is in a one- 
to-one correspondence with the quadratic UV divergence of the twist-4 operator. Taking the IR 
renormalon contibution as an estimate of real twist-4 matrix elements is therefore equivalent to 
the assumption that the latter are dominated by their UV divergent part. 

In practice, the all-orders calculation can be performed only in the large- Np limit, which 
corresponds to the QED case, where it produces an ultraviolet renormalon. To make the con- 



nection to QCD one uses the 'Naive Non Abelianization' (NNA) recipe |12| and substitutes for 
the QED one-loop /9-function the corresponding QCD expression. 

12 1 / 2 \ 

$, imit = N F -> — 11 - -N F . (12) 

The positive sign of (3$ in QCD produces an infrared-renormalon. While no satisfactory justi- 
fication for this approximation is known so far, it has proven to work quite well in low orders 
were comparison to the known exact coefficients in the MS scheme is possible ||13|| . Despite its 
phenomenological success in describing the shape of higher-twist corrections to various QCD 
observables, one should be aware of conceptual limitations of this approach. As the renormalon 
contribution is constructed from twist-2 parton distributions only, it really has no sensitivity 
to the intrinsic non-perturbative twist-4 nucleon structure. For example the ratio of the n-th 
moment M n of the twist-4 renormalon prediction to a structure function and the n-th moment of 
the twist-2 part of the structure functions itself is by construction the same for different hadrons: 



M. 



twist— A 



J^twist—2 



M, 



twist— 4 



hadronl n 



= , (13) 

hadron2 



which is definitely not necessarily the case in reality. 

As far as the scheme dependence of the renormalon model is concerned, its prediction for 
the magnitude of the mass scale of higher-twist corrections ~ A^e~ c is scheme invariant in the 
large- Np limit only. On the other hand, once one takes this scale to be a fit parameter, where 
the magnitude of this fit parameter depends on the process under consideration, this scheme 
dependence disappears. Even in this case, however, the mass scale still depends on the scheme 
used to determine perturbative corrections to the twist-2 part. Thus one can only say that the 
renormalon model is just a very economical way of estimating higher-twist contributions in the 
situation when their exact theory has not yet been constructed. In addition, in cases where the 
operator product expansion is not directly applicable, like in Drell Yan, the renormalon can still 
trace power corrections and give an approximation of their magnitude. 

In order to demonstrate how it works, let us consider the phenomenologically important case 
of twist-4 corrections to the Bjorken sum- rule: 

Cdxg{- n {xM 2 ) = \- + ^^ + 0{l/Q') 
Jo 6 g v Q z 



= -1.2601 ±0.0025 [14]. (14) 

The 1/Q 2 power corrections consist of target-mass corrections, a twist-3 and a twist-4 part ]Ej|. 
Only the twist-4 part can be traced by the renormalon ambiguity because this is the only oper- 
ator which can mix with the spin-1, twist-2 axial current operator. Table 1 contains results of 
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method 



value 



reference 



sum rule 



bag model 



chiral soliton model 



extraction from data 



renormalon 



fi = 0.003 ± 0.008 
Vhwnt-4 = -0.0094 

^ = 0.027 
Htwist-4 = 0.014 

fHwist-4 — —0.013 

Htwist-i = 0.012 ± 0.037 

H = ±0.017 



m 



ITf 
T 



m 
m 



Table 1: Higher-twist contributions to the Bjorken sum rule 



various approaches employed so far. A QCD-sum rule calculation is described in JT6|, Jj. Due 
to the high dimension of the corresponding operator the error arising from the sum rule calcu- 
lation is large, typically of the order of 50 percent. Another possibility is the MIT-Bag model, 
which however lacks explicit Lorentz invariance. As long as an unambiguous field-theoretical 
definition of twist-4 matrix elements does not exist, only the twist-3 part of the coefficient \x 
can be calculated on the lattice Recently the twist-4 matrix element has been estimated 
within the framework of the chiral soliton model [F8] , providing a semi classical description. 
The prediction of the renormalon model, which attributes the whole power correction to the 
twist-4 operator, is in an order-of-magnitude agreement with results of the other methods. The 
twist-4 contribution has also been estimated from the experimental data [BJ, yielding a value 
within a large error which is also in agreement with all the other methods as far as the order of 
magnitude is concerned, but predicting a positive sign in contrast to all calculations except the 
bag-model calculation. Note, that due to its nature as an ambiguity in the twist- 2 calculation 
the renormalon prediction is not capable to fix the absolute sign. 

Having the calculational efforts of sum rules and lattice QCD in mind, it might be an inter- 
esting idea to get a zero-order guess of the size of the twist-4 corrections by looking at the 
ambiguity in the expansion of the twist- 2 part. 



3 The renormalon contribution to F^^F^F^ and g\ 

Finally we will discuss the application of the renormalon model to the x dependence of DIS 
structure functions. Indeed the real success of the renormalon method in DIS - as noticed in 
|| 0| - is their capability to reproduce the x dependence of measured higher-twist contributions. 
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The absolute magnitude of the renormalon contribution can be either left as a fit parameter, 
taken from the NNA-calculation, or taken as an universal constant as done in the gluon scheme 



of 21 



We once more regard the nonsinglet part of the structure function g\ |22j. In Fig. [I] xg\ is 
plotted from the parton distribution set Gehrman/Stirling Gluon A |23|| plus/minus the renor- 
malon predicted twist-4 contribution, which is on the percent level. The dotted line gives the 
renormalon prediction again, amplified by a factor of 10. It becomes visible that the curve con- 
tains a zero at x ~ 0.7, a feature which is a very definite and clear prediction and should be 
confronted to experimental data, if a precision on the precent level was reached. 



As to Fl 0, we plot in Fig. |2| the appertaining renormalon prediction to the twist-4 opera- 
tor (short dashed line). The target mass corrections taken from are plotted by a dotted line. 
The sum of both is combined to the long dashed line and has to be compared to the experminen- 
tal fit by P5| , |2(| (solid line). The x dependence is quite well approximated by the renormalon 
approach, which predicts well the measured x dependence although the absolute magnitude is 
smaller than what is suggested by the experimental fit, but the x dependence is described in an 
acceptable way. 

As to the notation for the twist-4 corrections of the structure functions F2 and F3 we write 
(i = 2,3) 

F 4 (*,Q 2 ) = # T) (x,g 2 ) + ^(x) = Ff T )(x,g 2 )(l + ^ . (15) 

In the case of F 2 the good agreement between the x behavior of the nonsinglet renormalon con- 
tribution and the deuteron and proton twist-4 contribution (see Fig |3|) to F 2 has been noticed 
in . It is remarkable that the uncalculated pure singlet part seems to be small or proportional 
to the nonsinglet part. At least in the large-x range, where gluons contribute only a minor part, 
this is understandable. If we calculate the absolute value of the renormalon contributions for 
F 2 (solid line) it falls short by a factor 2 or 3, compared with what seems to be required by the 
data §. 

The same behavior is found for the structure function F3, where the twist-4 contribution is 
shown in Fig. f|. Again, if one is adjusting the absolute size of the renormalon contribution in 
the same way as it was done for F 2 , the result is in the large- x range in a good agreement with 
the data. Even more remarkable is the fact that the data indicate a change in sign in the large- a; 
region close to the one of the renormalon prediction. 



4 Conclusions 

The renormalon prescription provides a satisfactory description of the shape of the measured 
values of the twist-4 contributions at large x in the cases studied so far, leaving the absolute 
normalization as a fit parameter. Empirically, it falls short by a factor 2 - 3 as to the measured 
magnitude of the finite twist-4 contributions, when it is calculated in the MS scheme. 
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Figure 1: Experimental fit [|3[] for gi(x, 4 GeV 2 ) (full line). The estimate for the twist-4 correction 
is given explicitly multiplied by a factor of 10 (dotted line). The correction was added and subtracted 
from the experimental values (dashed lines) (A^ = 200 MeV, Q 2 = 4 GeV and N f = 4). 
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Figure 3: The solid line shows the renormalon model predictions for C p and C^, which are 
indistinguishable, LO parameterization taken from |27|. The dashed line shows the fit, as in 
Ref. [0. The filled and empty circles display the data for C p and Cd according to Ref. 0, 
respectively. 




Figure 4: Renormalon prediction for xh^{x) using the LO GRV |27j parametrization (solid line). 
The data points taken from Ref. [fZH] correspond to the NLO analysis, which provides the best 
fit of the higher-twist corrections. However, to be consitent within the renormalon approach 
the parton distribution sets used for the renormalon prediction have to be LO. The dashed line 
shows the prediction with the scale fi 2 adjusted to the description of the coefficients C p and Cd, 
as in H. The dot-dashed line shows the original prediction of f7|, obtained using the MRSA 
parametrization normalized at Q 2 = 10 GeV 2 . 
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